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HYDRODYNAMIC APPROACH TO CONSTRUCTING 
SOLUTIONS OF NONLINEAR SCHRODINGER EQUATION IN 
q i THE CRITICAL CASE 
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, Abstract. Proceeding from the hydrodynamic approach, we construct exact 

solutions to nonlinear Schrodinger equation with special properties. The so- 
lutions describe collapse, in finite time, and scattering, over infinite time, of 
wave packets. They generalize known blow-up solutions based on the "ground 
U» 1 , state " 

< 

-s. 

1. Preliminaries 

Consider the initial value problem for nonlinear Schrodinger equation (NLS) in 



> 
00 



C3 



l> 1 1 . 



W' t + A* + = 0, *(x,i):t;xI+^C, (1) 

*(0,x) = * (x)ei/ 1 (R"). (2) 

It is well known that the Cauchy problem (1),(2) has locally in time a solution 
of class C([0,T);H 1 (R n )), T < oo 0, g], El- The Cauchy problem is also locally 
well-posed in L 2 (R")|3], and this space is optimal 
■<^j- | Moreover, T = oo for a < ^, where dispersion dominates, and for a > ^ the 

solution may "blow up" in finite time under certain initial conditions, e.g.|S],[7]- 
"^h \ More exactly, there exist initial data and a positive constant T(5 , o( x )) < 00 , 

such that lim / |V*(x, t)\ 2 dx = 00. 

t-T(* (x))7 

The blow-up corresponds to self-trapping of beams in the laser propagation. 
A very good review with references can be found in |S] 
, For solutions to (1),(2) the following quantities are conserved: 

N[9] = ||*(x,t)|| ia(R «) = ||*o(x)|| L 2 (R „ ) = N (the probability), (7) 

P[\fr] = Im J *(x, i)V*(x, t)dx — P (the linear momentum), (II) 

and 

H[#] = HW^tJHi,^) - -^-||*( x ,i)||-+2 3(Rn) = H (the energy). (Ill) 
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Let u = — . Assume additionally a stronger decay of initial data at infinity, 
namely, |x|^'o( x ) G Li2(M. n ). Then the following important identity holds: 

M"{i) = 8H, (3) 

where M(t) = J |*(x, t)| 2 |x| 2 dx. This identity appears in many papers, but |U] 

seems to be the earliest one. It implies a very simple sufficient condition for the 
blow-up. It is easy to see that M (t) > for the solution of class C([0, T);H 1 (W 1 )), 
T < oo. At the same time (3) implies that 

M(t) = AHt 2 + M'(0)t + M (0). (4) 

The blow-up corresponds to vanishing of M(t). It signifies the concentration of the 
solution support in a set of zero measure. From elementary algebra arguments wc 
get that M(t) — > at a finite moment of time T if K < 0, where 

K = \QHM{Q) ~ (M'(0)) 2 . (5) 

This time T is positive if H < or H > 0, M'(0) < 0. 

In it was proved that for radial initial conditions ^o( x ) with H < the 
solution blows up in finite time in H 1 (M n ) with no integrability conditions on 
|*| 2 |x| 2 (see also ^T] in the context). 

2. Hydro dynamic interpretation 

Below we use the hydrodynamic approach due to Madelung, e.g. ^21- Namely, 
we represent the solution in trigonometric form, that is ^(t, x) = A(t, x) exp(i0(t, x)), 
where A(t,x) > and <fi(t,x) are real functions, the amplitude and the phase of 
the wave, respectively. 

Note that if A(t, x) is compactly supported, then x) is compactly supported, 
too. The support does not depend on the phase function <f>. 

Substituting this representation in (1) and taking the real and imaginary parts 
of the resulting equation, we obtain the following system 

A' t + 2{VA,V(/>)+AA(j) = 0, (6) 

A4>' t + A\V(f)\ 2 = AA + A a+1 . (7) 

Further, we multiply (6) by 2 A and apply the gradient operator to (7). Denote 
by p the probability density, A 2 = | ^> \ 2 , and by V double the gradient of the phase 
function, 2V</>. Thus, the final hydrodynamic form of (1) consists of two equations 

p' t + P divV+(V,Vp) = Q, (8) 

p(V t ' + (V,V)V) = 2(AV(AA + A a+1 )-VA(AA + A' J+1 j), A = p 1/2 . (9) 

The only difference from the traditional gas dynamics is the "exotic" pressure on 
the right-hand side of (9). This type of pressure changes the character of singularity 
completely, but allows to use the same methods as in the gas dynamics. 

The data 

p(0,x)=po(x), V(0,x) = V (x) 
complete the statement of the Cauchy problem for (8), (9). 
It suffices to demand that 

^ /2 (x) G H l (B»), pI /2 (x)V (x) G L 2 (R n ), 
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to ensure that the corresponding initial function \Po(x) belongs to the class H 1 (R n ). 
Note also that p vanishes as |x| — * oo, however, the same is not necessary for V. 
The conservation laws (I— III) in the new terms are as follows: 

Jpdx = N (the mass), (I 7 ) 
J pVdx — P (the linear momentum) (-^0 



and 



^ + ^--r^ /2+1 ) dx = H (the energy). {Ill') 
4 4p a + 2 



One can see that the "kinetic energy" component in H is the same as in gas 
dynamics (up to the multiplier). 

Note that in gas dynamics terms we have 

M(t) = J p(x,i)|x| 2 dx, 

by virtue of (8) 

M'(t) = 2 /"(x,V(x,f))p(x,f)dx ) 



here x is the radius-vector of point in the space. 
Let us introduce one more functional: 



Q A {t) = j '(x,A)p(x,i)dx, 



where A is a constant vector from M. n . From (8) and the conservation of linear 
momentum (II') we have Q'^(t) — (P, A) := Pa = const, and 

Qa(*)=Pa* + Qa(0). (10) 

From the Holder inequality we have also 

Q\{t) < \A\ 2 NM(t). 

In the domains where the amplitude A > 0, instead of (9) we consider the 
equivalent equation 

AA 

Vl + {V,V)V = 2V(— +A°). 

Now we use the following idea, recently applied to construct solutions to the gas 
dynamics equations fsee |13| . |14j . |15|. [lbp. Namely, let the "velocity field" have 
the form 

V = a(t)x, (11) 

where alt) is a time-dependent function. Thus, the phase function can be restored 
as 

0(t,x) =o(t)i^-+7(t), 
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with an unknown function 7(i). Then, from the (linear in p) equation (8) we find 
the density as 

p(t,x) = exp(— n j a(T)dT)p I xexp(— J a(r)rfr) J , 



or 



A(t,x) = exp(-^ J a(r)dr)A I xexp(- ^ a(r)dr) J , 



with A (x) = |* (x)|. 

Note that in the critical case, <r = -, 



AA(t,x) 



t 

+ ^(t,x)=exp(-2 J a (r)dr)(^|)+^(0), (12) 



where £ = xexp(— J" a(r)dr). 

o 

Further, we have from (7) and (12) 



t 

(a'(i) + a\t))^ + j'(t) = cxp(-2 J o(r)dr) + A„ 



(0 



or 



t t 
(fl '(t) +fl2 (t))M! + y W exp ( -2/ a( ^).exp(-4/ fl( ^)( 



AA (g) 

4>(0 



The variables t and £ can be separated if 



r 

/(t) =7oexp(-2 y o(t)), 70 el 1 . 



In this case 



(a'(t) + a 2 (t))|e| 2 = 4exp(-4 1 o(r)dr) (^|^ + ^o(0 - 7o 

o 

It follows from (14) that 

a'(t) + a 2 (t) = 4fccxp(-4 / a(r)dr), 



AAq(x) + ^ +1 (x) = (fc|x| 2 + 7 o)A,(x), 

where fc is a constant. 

Thus, we seek a special solution to (1) in the form 

t r t 



(13) 
(14) 

(15) 
(16) 



5'(t, x) = exp(— — y a(r)dr)Ao xcxp(— J a(r)dr) 
o L o 

cxp ^ia(t)^J-^ CX P ^*7o y cxp(-2 y a(Ti)dri)dr^ exp(i#), (17) 
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with flel 1 . 

Because the origin is not a particular point, without loss in generality we may 
consider the velocity field V = a(t)(x — x ), where x is an arbitrary fixed point. 
Then the solution takes the form 

t r t 

*(t,x) = cxp(-^ J a(T)dr)Ao (x - x ) exp(- J a(r)dr) 



cxp ^ia(t)— — ^ X °^ ^ cxp ^ijo J exp(— 2 J a(ri)dri)rfr^ cxp(z#). (17') 

The decay properties of the solution as |x| — ► oo depend on ^4o( x )- The physical 
sense requires that the solution should be of the class L 2 (R"). If we wish to consider 
solutions from the space C([0, T); H 1 (R n )) 7 T < oo, natural for the existence and 
uniqueness to the Cauchy problem (1), (2), we have to choose the initial data 

/ |x — x | 2 \ 
*o(x) = A (x - x ) exp I ia exp(i9), 



where Ao(x) is a non-negative solution to (14) belonging to _ff 1 (IR™) and |x|Ao(x) e 
L,2(M. n ). Note that with the function Aq of this class we have also conservation laws 
(I - III) for solutions (17) and (17'). 

3. Time evolution 

Let us investigate the qualitative behavior of a(t) governed by (15). It is easy 
to see that in the case k < for any initial datum a(0) = ao there exists a moment 
T* such that a(t) — > — oo, t — > T*. Really, as a'(t) < e < 0, then under any a 
there exists a moment T\ > when a(Ti) < 0. Further, the comparison theorem 
shows that a(t) < a(t), where a(t) is a solution to the Cauchy problem a{t)' = 
-5 2 (t), a(Ti) = a(Ti), i > Ti. This means that if fc < 0, the solution * to (1) of 
the form (17), with the amplitude Aq satisfying (16), localizes at the origin within 
a finite interval of time, provided N < oo. 

We can also analyze the differential corollary of (15) 

a"{t) +6a(t)a'(t) + Aa 3 {t) = 0, 

with the initial data a(0) = a , a'(0) = — Oq + 4k. It gives us, in particular, that in 
the case k > the solution vanishes at infinity as 0(i -1 ). 

If k = 0, any nontrivial solution to (15) blows up at a finite time T = —a^ . 
The time is positive if a < 0. Moreover, if |x|\I>o( x ) <= L 2 (M. n ), we can find a(t) 
explicitly for any k. Observe that, if V = a(t)x, then 

M'(t) = 2a(t)M(t), (18) 
therefore a(t) — ^g^y- The explicit form of M(t) is known, sec (4). Thus, 

%Ht + M'{0) 

° W ~ 2{AHt 2 + M'(0)t + M(0))' ( ' 

If K = 16HM(0) - (M'(0)) 2 < (sec (5)), then M(t) tends to zero (a(t) goes to 
infinity, respectively) within a finite interval of time. Moreover, this interval can 
be readily calculated. If K > 0, then a(t) ~ t — > oo. 
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Now taking into account (15), (18), (19) we can compute 

WHMjfi) - (M'(0)) 2 f2m 
16M 2 (0) 16M 2 (0). 1 

There are situations where we can express a(t) through Q\(t)- Namely, for V = 
a(t)x we get P A = a(t)Q A (t). Therefore, if there is A e 1" such that Q A (t) ^ 
(Pa 7^ 0), then taking into account (10) we have 

1 1 



a(t) 



i + (?A(0)P A 1 t + a Q 



l • 



From (18) we obtain now M(t) = M(Q)al(t + a^ 1 ) 2 . Comparing the result with (4) 
we can see that AH = Af(0)a 2 (0) and K = k = 0. 
Note that if A is radial, then Q\(t) = for any A. 

4. Evolution of wave packets 
Summarizing the above results, we can formulate the following theorem: 

Theorem 1. Suppose that (16) has at least one nonnegative solution A (x). Then 
equation (1) has a special solution given by the explicit formula (17)((1T)), with 
the function a(t) governed by equation (15). 

IfxA (x) G L 2 (R n ), then the formula (IT) can be written as 



, , /M(0)\ 4 „ 



\M{t)) 



(x- x ) 



exp ^^^|x-x | 2 ^ exp (i 7o M(0) J M-\r)dr^j exp(i0), (21) 

with the quadratic function M{t) having the explicit form (4). 

The behavior of the solution depends on the sign of the constant k (see (20)). 
If k > 0, then the solution decays. Namely, 

max |*(t,x)| =A+ ( MM] 4 t -»/2 A+ = maxA (x). 
xei" y AH ) r™ 

If k < 0, then the solution blows up at the point x at a finite moment of time 
T. This time is positive in the following cases: 
(i) ifH = 0, M'(0) < 0, then T = 

-m'(o)-tJ\k~\ 



(ii) ifH>0, M'(0) < 0, then T - - 7i7 

(in) ifH<0, thenT= Z^pM . 
Moreover, for k = 0, 



for k < 0, 



max|*(t,x)| =A+ ( ^BA * {T-t)- n '\ 
xei" \ J\K\) 



HYDRODYNAMIC APPROACH TO CONSTRUCTING SOLUTIONS OF NLS 



7 



Remark 1. If there exists Ael° such that Q\{t) ^ (Pa ^ 0), then only the 
situation with k = may be realized. 

Remark 2. For the solution of form (19) we have 

l|V*(x ) t)||| 2(R „ ) =exp(-2 J a(*)dt)||VA)(x)||| 2(R n ) + 
ia 2 (i)exp(2 J a(t)dt)\\xA {x)\\] 



^||VAo(x)|| i2(R „ ) + T ^--0(T-t) , 

where A = 1 (the lower estimate for the blow-up order, see f° r k < 0, and 
A = 2 for k = 0. 

5. Comparison with previous results 

In the theory of NLS for the critical case a — — the crucial role is played by the 
so-called ground state, i.e. the positive radially-decreasing solution to the elliptic 
problem 

Au+\u\ a u-u = 0, ueH l (R n ). (22) 

It is known ^H] , that the solution with such properties is unique and exists at least 
for n = 1,2,3. The solution belongs to C 2 (R"), and \D a u\ < Cexp(-S\x\), where 
C, S are positive constants. We denote it by R(x). 
It is known that H[R(x)} = 0P|. 

It was proved ^Hl that if ||5'o( x )||l 2 (R") < II-R( x )Il 2 (R")> then the solution to the 
problem (1), (2) is global in time. If the solution blows up, then ||\1/o(I x )||l2(k™) > 
||i2(x)|| £ 2 (M»)- In the case ||^'o( x )||l 2 (R") = II-R( x )IIl 2 (r») the solution either blows 
up or not. 

If it does blow up, it necessarily has the following special form based on the 
ground state |2Till5T] : 

n / t\\ ( , |x — Xnl 2 A / UJ \ T „ / Uj(x — Xn) 

*(x,t) = e X p^)exp^(— + L_^)j (__) *(U__£l_ Xl 

(23) 

with certain x el",Xi el",SeM,ueR + ,TeR. However, (23) coincides with 
(21) for K = k = and 70 = 1. Indeed, in this case R(x) is a solution to (14), 

M (t) = H(t + ^l) 2 =H(t (-^P) ") 2 , lo 2 = MM, T = uj= Here we 

choose M'(0) < to guarantee the positivity of T. Note that H > for solutions 
of the form (23). 

It can be readily demonstrated that 

Af(0) 2 2 {M'{t)f 

^(0) , , (M'(t)f 

Thus, for A = R, we get H = ^M(t) > °- 

It follows from the above results that if Aq(x) is a solution to (16) with k < 
from H 1 (M n ), then 

||Ao( x )|| L 2 (R „) > ||i2(x)|| £a(H n), 
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because in this case the solution (17) blows up and it is not of the form (23). 

The profile of our solution in the case k < is different from that of the solution 
(23) (corresponding to k = 0), and so is the rate of blow-up (see Remark 2 of 
Section 4.) 

Note that in the case k < the solution to equation (16), considered in the 
space R™, oscillates as | as j — > oo. For linearized equation (n = 1) we can even get 
the explicit solution: 



u = | a; |-V^C 1 W 1 (^|=, -, -W2kx<) + C 2 W 2 (^=, -, 2»V2fc^)), 

where W\^W 2 are the Whittaker functions, C^C^ are constants. This highly os- 
cillating function does not belong to L 2 (R ra ). So we cannot hope that the solution 
to nonlinear perturbed problem (16) is positive and belongs to L 2 (M. n )(H 1 (M n )). 

However, we can consider (in higher dimensions too) the solution to (16) given 
in £1 C R™. To be exact, now we deal with the Dirichlet problem for (16) with zero 
boundary conditions. Regarding this Dirichlet problem, for example, there is the 
following result due to [2"2"j . 

Theorem 2 ( 22 ). Let be a bounded, smooth domain in R™,n > 2, and g : 
R + x f2 — * R + a locally Lipshitzian map. Consider the elliptical boundary problem 

Am + u p + eg(x,u) — 0, u\qq — 0. (*) 

If 1 < p < (p > 1 for n = 1, 2), there exists eo > such that, for < e < €q, 
(*) has a solution u — u e , which is positive on Q. 

_Ul |2_ 

In our situation g(x, u) = — — — —u satisfies to the theorem condition for k < 
in a ball from R" , for 70 < this condition holds for all l n . 

The solution to the Dirichlet problem is classical, that is it belongs to C 2 (f2) D 
C(f2), therefore it can be extended to all over the space R" at least as solution from 
L 2 (R"). 

6. Further generalization 
Let us consider the following velocity field: 

V = a(t)x + b(t)A, (24) 
with a constant vector A. Then the phase function takes the form 

0(t,x) = a(t)^- + h(t)(A 7 x)+ 1 (t). 

Proceeding in the spirit of Section 2, we obtain 

t / t t 

A{t, x) = exp(— ^ J a(r)dT)A I xexp(— J a(r)dr) — A J 6(r)exp(— J 

V 

t t 

Denote £ = xexp(— J a(r)dr) — A f Q 6(r) exp(— f r a(Ti)dri)dr. From (7) and 


(12) we get 

[a'(t) + a 2 (t)] M!+ 

(a'(t)+a 2 (t)) f b(r)exp(- [ a^dr^dr + (b' (t) + a(t)b(t)) exp(- / a(r)dr)' ^ A 
Jo Jo Jo 



a(ri)dri)ciT 
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(a'(t) +a 



Jo 



6(r)exp(- / a(Ti)dTi)dTy + 



2(b'(t) + a(i)6(i))exp(- f a(r)dr) / b(r)exp(- f a(T 1 )dT 1 )dr+ 
Jo Jo Jo 



6 2 (i)exp( 



-2 [ a 
Jo 



(r)dr) 



t t 
- 7 '(t)exp(-2 1 a(r)dr) = cxp(-4 1 a(r)dr) + ^o(o) • 

o o 
There are two possibilities for the separation of variables. 

I. The functions a(t),j(t) and Ao(xi), xi = x + A satisfy equations (15), (13) 
and (16), respectively, b(t) = a(t). So we return to the formula (17') considered 
above. 

II. The functions a(t),b(t),j(t) and A (x) satisfy the following equations: 

a'(t) + a 2 {t) = 0, 



b'(t) + a(t)b(t) = 2fci cxp(-3 / o(r)dr), 

Jo 

Y(t) =7oexp(-2 / a(r)dr)- 
Jo 

/»t 

4fciexp(-2 / a(r)dr) + / 6(r)exp(- / o(ri)dri)dr + 6 2 (i) 
Jo Jo Jo 



(25) 
(26) 

(27) 
(28) 



AAo(x) + ^ +1 (x) - (fci(A,x) + 7 o)Ao(x), 

where 70 and k\ are constants. 

The functions a(t),b(t),j(t) can be found explicitly. 

Note that if k\ = 0, 70 = 1, fr(i) = 0, we obtain again a solution that can be 
represented by formula (23). 

The simplest situation is a(t) = 0, b(t) = 2k x t + b , j(t) = -j^^fat + b ) 3 + 
7o* + 71, 7o, 71 € W\ The corresponding solution has the form 

tf (t, x) =Ao[x- At(kit + b )} cxp (^(fat + y )(A, x) + 7 (i) 

with the function ^4o( x ) satisfying (28). For fci = 0, 70 = 1, 60 = we get the 
solitary wave solution R(t, x) exp(z(i + 71)). 
In the general case 

a(t) = 7T i =Tl 6(t)= b ° a ° 1+ y -2^- 1 



7(0 



t + a 



To 



/,vA| 2 ? b(n) 



(t + O 2 ' 



(t + a 



ti + a 



T dn - -|A|V(r) 



dr. 



The corresponding solution has the form 



*(x, t) = exp(*7(t)) exp ( i(b(t)^l + 



4(t + a " ) 
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ll^T A /-^-r'/r- (2!)) 
t + a J t + a Q 

with the function A (x) satisfying (28). The solution (29) has an interesting feature 
in the case when k\ ^ 0, An(x) e L,2(M. n ) and ao < 0. At a finite time, T = — ajj" 
the mass concentrates at a point which escapes to infinity. 

If we suppose that A (x)x S ^QR™), we can express the solution through the 
functionals M(i), Q A (t) and P. Indeed, for the velocity field (24) we have 

M'(t) = 2a(t)M(t) + b(t)Q A (t), (30) 

Pa =a(t)g A (t) + fo(^|A| 2 . (31) 
Further, from (3), (30), (31) we obtain 

a'(t) = a(t) 2pAQA(t) ~ M '® N \ A \ 2 P A - 8HN\A\* 



M(t)N\A\*-Ql(t) M(t)N\A\*-Ql(ty 
Then we use (4) and (10) to get 

*®-^im + m (32) 

where F(t) = At 2 + Bt + C,A = 8HN\A\ 2 — P|, B = M'(0)HN\A\ 2 - 2P A Q A (0), 
C = M(0)7V|A| 2 - Q A (0). Note that by virtue of the Holder inequality C > 0. 
Further, we find from (32) and (31) 

M Ct , u .P A (A-C)t 2 + (P A (B-ao)-CQ A (0))t+(P A C-aoQA(0)) 
a(t) -fW +fl °' Ht) - F(t)N\M 2 ' 
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